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Abstract 

We have performed exhaustive multicanonical Monte Carlo simulations of three 12-residue syn- 
thetic peptides in order to investigate the thermodynamic and structural properties as well as the 
characteristic helix-coil transitions. In these studies, we employ a realistic model where the inter- 
actions between all atoms are taken into account. Effects of solvation are also simulated by using 
an implicit solvent model. 
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I. INTRODUCTION 



It is well known that three-dimensional (3D) conformations of peptides and proteins 
play an important role due to their biological activities. Therefore, many theoretical and 
experimental studies focus on the determination of the 3D structure of these molecules. In 
a newly growing field of research, synthetic peptides are investigated for use in hybrid nano- 
devices, depending on their self-assembly properties . In these studies, it is also shown 
that the binding of peptides on metal and semiconductor surfaces depends on the types of 
amino acids [3| and on the sequences of the residues in the peptide chain . These 

experiments reveal many different interesting and important problems, which are related to 
general aspects of the question why and how proteins fold. This regards, for example, the 
character of the adsorption process, i.e., whether the peptides simply dock to the substrate 
without noticeable structural changes or whether they perform conformational transitions 
before binding. A related question is how secondary structures of peptide folds in the bulk 
influence the binding behavior to substrates. In helical structures, for example, side chains 
are radially directed and - due to the helical symmetry - residues with a certain distance in 
the sequence arrange linearly. This could have consequences for docking to a regular crystal 
surface, where the atoms are also arranged linearly along the main axes. This means that two 
peptides with the same content of residues, but different sequences, could exhibit completely 
different binding properties. This behavior was, in fact, observed in a recent experimental 
adsorption study of peptides in the vicinity of semiconductor substrates , although other 
explanations for this kind of specificity are also conceivable. For these reasons, it is likely 
that the binding specificity also depends on the thermodynamic and structural properties of 
peptides in solvent, as it has already turned out in investigations of minimalistic models j?, ^. 

In this paper, we focus on three synthetic peptides, AQNPSDNNTHTH, AQNPS- 
DNNTATA. and TNHDHSNAPTNQ H. whose b,„di„g properties were invesfgated .n 
recent experimental studies [5|, \^. The second sequence is a mutated version of the first 
one, where the histidine residues (H) of the first chain are replaced by alanines (A). The 
third sequence is a randomly permuted sequence of the first chain. It was shown in the ex- 
periments that the first sequence has a strong affinity to adsorb to gallium arsenide (GaAs), 
whereas the binding to silicon (Si) is very weak. Exchanging the histidines by alanines 
improves the binding properties to Si, while the adsorption strength to GaAs is noticeably 
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reduced. For the randomly permuted sequence, the binding strength to GaAs is left widely 
unaffected, while binding to Si is as strong as to GaAs. 

Employing multicanonical (MUCA) Monte Carlo sampling we analyze single- 

molecule folds in the bulk and illuminate the thermodynamic and structural properties of 
these peptides. Generalised-ensemble methods applied to all-atom descriptions of proteins 
have been very successful in the past, e.g., in revealing the statistical mechanics in the 



folding process of small proteins 
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14l |. For sequences with more than 20 residues. 



studies of thermodynamics and kinetics employing realistic physical models are computa- 
tionally extremely demanding. For such systems, reduced all-atom models [l^ or models at 



a higher coarse-grained level [l6|, ll7| could be, depending on the particular question, much 



more promising. Coarse-grained lattice and off-lattice models allow for systematic thermo- 
dynamic studies and, at least partly, sequence analyses of heteropolymers with up to 100 
monomers [l8, 19, 20]. 



Firstly, we have simulated all three molecules in vacuum. Then, in order to see the 
effects of solvation, we have also performed extensive simulations of the commonly used 
surface-accessible area solvent model with OONS atomic solvation parameter set 2l|]. The 
preferential properties of this parameter set compared to others were reported in previous 



works 



22 



23|. 



The rest of the paper is organized as follows. After the description of the peptide model 
and the computational methods used in this study in Sect. [Tll we discuss in Sect. IIIII ex- 
emplified effects of solvation compared with results obtained in the vacuum simulations. In 
Sect. IIVI we discuss in detail the helix-coil transitions for the three peptides in solvent by 
means of fluctuations of several energetic and structural quantities. Section |V] addresses 
the folding channels in the free-energy landscape. The paper concludes with a summary in 
Sect. IVIl 

II. PEPTIDE MODEL AND SIMULATION METHOD 



A. The Peptide Model 

In our simulations a peptide is modeled with all of its atoms. Each atom i, located at 
the position r^, carries a partial charge qi. Covalent bonds between atoms, according to 
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the chemical structure of the amino acids, are considered rigid, i.e., bond lengths are kept 
constant, as well as bond angles between covalent bonds and certain rigid torsion angles. 
Distances between nonbonded atoms i and j are defined as Vij = |rj — rj\ and measured 
in A in the following. The set of degrees of freedom covers all dihedral torsion angles 
$, = {^a} of ath residue's backbone (0^, V'a, ^a) and side chain {x = Xa\xa\ ■ ■ ■)■ The 
model incorporates electrostatic Coulomb interactions between the partial atomic charges 
(all energies in kcal/mol), 

^c(«- 332i:-^. (1) 



effective atomic dipole-dipole interaction modeled via Lennard- Jones potentials 

4 ■ . R. • 



0-H and N-H hydrogen-bond formation. 



and considers dihedral torsional barriers (if any): 

EU^) = Y,Ui{l±cos{niCi)). (4) 
I 

The total energy of a conformation, whose structure is completely defined by the set of 
dihedral angles ^, is 

Eo{$) = Ec{^) + Elj(^) + ^hb(^) + ^tor(^). (5) 



The parameters qi, Aij, Bij,Cij, Dij,Ui, and ni are taken from the EC 



]PP/3 (Empirical 

Conformational Energies for Proteins and Polypeptides) force field 25|], one of the most 
commonly used all-atom force fields. In all simulations the dielectric constant was set to 
e = 2, which is the vacuum value. Proline's is considered rigid at —68.8°. We always 
used the trans down-puckering conformation of the proline ring. For the implicit-solvent 
simulations, the model is extended by the solvation-energy contribution, which is given 



by 



26| 



^solv(0 = E^^^*(^)' (6) 



where Ai is the solvent-accessible surface area of the ith atom for a given conformation and 
CTj is the solvation parameter for the ith atom. The values for depend on the type of the 



TABLE I: Peptide sequences studied in this work. 



SI 
S2 
S3 



AQNPSDNNTHTH 
AQNPSDNNTATA 
TNHDHSNAPTNQ 



ith atom and are parameterized according to the suggestions given in Ref. 2l(|. The total 
potential energy of the molecule then reads 



Etot(0 = ^0(^)+^solv(0- 



(7) 



The described peptide model and the ECEPP/3 parameterization is implemented in the 
software package SMMP which we used for our study. 

In Table [U we have listed the three sequences investigated with this model. 



B. Multicanonical Sampling 

Multicanonical sampling 1^, UJ, |28j] is a generalized-ensemble method, in which con- 
formations are ideally sampled according to a fiat energy distribution Pmuca(-£') = const., 
i.e., the Markovian dynamics of the algorithm corresponds to a random walk in en- 
ergy space. The desired canonical distribution at a certain temperature T is given by 
Pca.n{E,T) ~ n{E) exp {—E/ RT), where n{E) is the density of states and the gas con- 
stant takes the value R ~ 1.99 x 10^^ kcal/Kmol in the units used in this paper. Since 
canonical and multicanonical energy distributions are trivially related via pcan{E,T) ~ 
W~^^^{E) exp {—E/RT) p^^ca{E), the main task is a precise determination of the multi- 
canonical weights Wuinca.{E) ~ n^^{E). 

The implementation of MUCA is not straightforward as the density of states n{E) is 
unknown a priori. Therefore, the weights Wj^^ca{E) have to be determined in the first 
stage of the simulation process by an iterative procedure until the multicanonical histogram 
H{E) ^ const, in the desired energy interval. An efficient, error-weighted estimation method 
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. We note that the ef- 



for the multicanonical weights is described in detail in Refs. 
ficiency of the determination of the multicanonical weights usually depends on the choice 
of the simulation temperature, which was in the present study Tsim = 1 OOOK. The reason 
is that, since the "fiat" energy histogram covers a larger region in subsequent recursions. 



energetic states are hit for the first time, where the multicanonical weights are still undeter- 
mined. Because the ratio of the weights controls the acceptance of a conformational update, 
the dynamics of the recursion part of the algorithm is noticeably influenced. This behavior 
can be "smoothed" by a careful choice of the simulation temperature. 

Eventually, in the second stage of the multicanonical simulation, a long production run is 
performed based on fixed multicanonical weights. Since the weights for all energetic states 
in the desired energy range have already been determined, the choice of a certain simulation 
temperature is unnecessary. 

In our concrete implementation, we first carried out Metropolis simulations at relatively 
high simulation temperatures and MUCA test runs which enabled us to determine the 
required energy range. This interval was then divided into bins of 1 kcal/mol. At each update 
step, a trial conformation was obtained by changing a dihedral angle C,i —>■ C,i within the range 
[—180°, 180°], which was accepted according to the transition probability uj{$, = 
mm[exp {S {E {$,))- S{E{^'))) ,1], where S{E{$,)) = -\ogW^nUE{^)) can be identified 
with the microcanonical entropy. The dihedral angles were always visited in a predefined, 
sequential order, i.e., a sweep is a cycle of N Monte Carlo steps {N = total number of 
dihedral angles). 

The weights were built in 200 recursions during a long single simulation, where the 
multicanonical parameters were iterated every 10 000 sweeps. Then, we performed a full 
simulation of two million sweeps with fixed weights, which covers the temperature region up 
to Tmax = 1 000 K reliably. In Fig. [1], the density of states and the multicanonical histogram 
for the first sequence considered, AQNPSDNNTHTH (SI), is shown. As seen from this 
figure, the multicanonical histogram is indeed "fiat" which is a necessary condition for 
the multicanonical technique to be reliably working. Statistical expectation values for any 
thermodynamic quantity A and all temperatures can finally be calculated from the time 
series recorded during the multicanonical production run: 

where labels the conformation at "time" t and (3 = 1/ RT is the inverse thermal energy. 
The derivative of the quantity (A) with respect to the thermal energy is given by 

d{A) 1 



d{RT) {RTf 
6 



{{EA)-{E){A)). (9) 
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FIG. 1: Natural logarithm of the density of states n(E) and multicanonical histogram from the 
simulation of the sequence AQNPSDNNTHTH (SI) in solvent. 




FIG. 2: (Color online) Low-energy reference conformations of the peptide SI in vacuum (left) and 
solvent (right) for the calculation of the overlap parameter (jlip . 



Expressions like this are typically used to discuss the influence of thermal fluctuations on 
A. 



III. SOLVATION EFFECTS 



In this section, we compare the folding behavior of the exemplifled peptide SI in vac- 
uum and solvent, respectively. Changes of energy fluctuations, i.e., the widths of the en- 
ergy distributions, signalize typically a crossover or transition between signiflcantly different 
macrostates ("phases") of the system considered. For polymers or peptides, the crossover 
between such macrostates is accompanied by a cooperative conformational transition. Thus 
it is reasonable to compare the behavior of the peptide in vacuum and solvent with regard 
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FIG. 3: Average angular overlap Q between ensemble conformations and the reference conforma- 
tions in Fig. [2] for the peptide SI as a function of temperature in vacuum and solvent. 



to energetic fluctuations, the specific heat (in units of R) 

1 



V 



29 



30 



{RTf 

and with respect to the angular overlap parameter 
("order" parameter), which is suitably defined as 

1 ^ 



(10) 



3l| as a structural quantity 



(11) 



i=l 



where d{a,a') = min(|Q; — q;'|,360° — \a — a'\). In this expression, the dihedral angles 
of the actual conformation are compared with the corresponding torsion angles of a 
suitable reference conformation, ^'^'^^ . Since the overlap parameter ( ITTi) is a measure for 
the similarity of any conformation and the reference conformation, it can be considered 
as a system state parameter: Q = 1 only if the considered conformation is identical with 
the reference conformation, which is here chosen to be a typical representative of the (low- 
energy) helical phase. In Fig. O the helical low-energy reference conformations of sequence 
SI in vacuum (left) and solvent (right) are shown. The main differences regard the side- 
chains and the tail, whereas the helical part is hardly infiuenced by solvation effects. The 
overlap Q between these two conformations comparing all dihedral angles is 0.595, while 
considering only the backbone dihedrals, Q ~ 0.690. In the latter case, the comparatively 
still small coincidence is mainly due to the non-helical tails, which are highly flexible, i.e.. 
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FIG. 4: Specific heat vs. temperature for the peptide SI in vacuum and solvent. 
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FIG. 5: Specific heat vs. temperature for the three peptides in solvent. 



there are no stabilizing hydrogen bonds with the backbone. Proline breaks the helix in both 
cases. 

In Fig. [3] we have plotted the average overlap parameter (Q) for this exemplified sequence 
in vacuum and solvent. For low temperatures, in both cases most of the conformations in the 
ensemble have similarities with the reference conformations, i.e., the ensemble is dominated 
by helical conformations. For the peptide in solvent, the average overlap parameter decreases 
rapidly at about 440 K, the conformations in the high-temperature phase are random coils. 
The situation is comparable for the peptide in vacuum, with the noticeable differences that 
(Q) decreases much slower at a transition temperature near 530 K. This is confirmed by 
considering the energetic fluctuations of the system, i.e., the specific heat as shown in Fig. HI 
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FIG. 6: Average number of helical residues (a) and derivative (b) vs. temperature for the systems 
in solvent. 

In the vicinity of the peak temperatures, the peptides exhibit conformational activity. The 
helix-coil transition peak is stronger and sharper for the peptide in solvent, the transition 
temperatures are close to the above-mentioned values. 

Summarizing, the main effect of the solvent is the strengthening of the helix-coil transition 
which is also present in the vacuum case. Furthermore, the transition temperature is shifted 
by about 100 K towards lower temperatures. These results are as expected, since it is known 
that solvent stabilizes secondary structures and therefore the barrier to resolve the helix is 
higher than in the vacuum case and the relaxation of the fluctuations of the peptide-solvent 
coupling degrees of freedom leads to a lower transition temperature. These differences have 



also been observed in studies using other parameter sets 



It should be noted that also in the OONS implicit-solvent model the transition tem- 
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perature is probably still strongly overestimated as is already known from studies of other 
helical peptides |33|. One of the reasons is the choice of a temperature-independent solvent- 
peptide coupling strength and the "smeared" , nonlocal and static polar environment without 
intrinsic fluid properties. 



IV. HELIX-COIL TRANSITIONS OF THE PEPTIDES IN SOLVENT 



In the following we discuss the thermodynamic properties of the three synthetic sequences 
employing the ECEPP/3 force field with OONS implicit-solvent parameter set. Considering 



several quantities we find strong indications for helix-coil transitions in all three cases. He^ 
coil transitions in peptides and nucleic acids were first addressed by Zimm and Bragg 
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and have been studied extensively 
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As a first indication for the conformational transitions, we find strong peaks in the specific 
heat as shown for the three peptides in solvent in Fig. [5l The peaks are located near 
440, 420, and 410 K for the wild-type SI, mutant S2, and randomly permuted sequence 
S3, respectively. The character of the conformational transition is identified by measuring 
the temperature-dependence of a-helicity and /5-sheetness. The helicity is a natural order 
parameter for the identification of helix-coil transitions in peptides. A residue is defined to 
be in a-helical state, if its backbone dihedral angles and are in the range (—70 ± 20°) 
and (—37 ±20°), respectively. We have shown the changes of this quantity and its derivative 
versus temperature in Figs. [6]^a) and (b), respectively. 

In the calculation of the helicity, end residues are not taken into account, because end 
residues are very fiexible and do not conform to any definite state. There is also a proline 
residue in all chains which is known as helix-breaker, because it lacks a primary amine group 
and due to the peptide bond to the preceding amino acid, there is no H-atom allowing for 
the formation of a hydrogen bond that could stabilize an a-helix or a /3-sheet structure. 
Furthermore, the rigid proline side chain typically forces for steric reasons the angle of the 
preceding amino acid to take non-helical values. Hence, the maximum number of residues in 
a helical segment is eight for the three sequences under consideration. Rapid decreases of the 
average helicities (as shown in Fig. [6](a)) are noticed for all three peptides. The transition 
temperatures lie in the same temperature region as the peaks of the specific heat, as can be 
read off from the fluctuations around the average helicity as plotted in Fig. [6](b). 
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FIG. 7: Average number of /3-sheet structures (a) and derivative (b) vs. temperature for the 
peptides in solvent. 

In the force-field parameterization used, /? sheets can be excluded in the low-temperature 
region. This is confirmed by the plots in Figs. [Tl^a) and (b), where the sheetness (^/j), which 
is the average number of Ramachandran angles in a /5-sheet state [i.e., G (—150 ± 30°) 
and G (150 ± 30°)], and the fluctuations of this quantity are shown. No noticeable 
/3-sheet structure is identified in the low-temperature region because the whole ensemble 
consists of strongly helical conformations. The average sheetness increases slightly above 
the conformational-transition temperature, but this signal is relatively weak and the high- 
temperature ensemble is expected to be dominated by random-coil structures. 

Another interesting quantity is the probability of the individual residues to become "he- 
lical". In Figs. [H](a)-(c), we have plotted the color-coded profiles of the probabilities of each 
residue to be in a helical state. Since it is not intuitive that a single residue can form an 
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a-helix motif, although its dihedral angles are in the range of a-helical region of the Ra- 
machandran map, we define a residue to be helical only if it is part of a helical segment with 
at least three successive helical pairs of Ramachandran angles. This allows a clearer view 
on the helix formation. For the three considered peptides in implicit solvent, the helix-coil 
transition is a sharp one-step process in the temperature region between 400 and 500 K. The 
single hehcal segment of the peptides SI and S2 is formed by the residues 5 to 12 (counting 
from the N terminus). Proline at position 4 in the sequence breaks the helix and the 1-4 
residual tail is coil-like. Surprisingly, the randomly permuted sequence S3, where proline is 
located at position 9 (which is, unfortunately, also at the fourth position - counted from the 
C terminus), exhibits in addition to the strong 1-8 helix a second helical segment between 
residues 10 and 12. Although this signal is weak, a non-negligible subset of conformations in 
the low-temperature ensemble contains two independent helical segments, broken by proline. 
The transition temperature for the 10-12 helix is slightly smaller than for the main segment 
and lies below 400 K. Note that the formation of the second helix for sequence S3 can also 
be observed as "shoulder" in the corresponding specific heat in Fig. [5l 

Hydrogen bonds are mainly responsible for the formation and stability of secondary 
structures such as a- helices. In Figs. IHt^a) and (b), the average numbers of hydrogen bonds 
(nn-bond) and their fluctuations, respectively, are shown for the three sequences. In the he- 
lical phase, conformations typically possess approximately 4-6 hydrogen bonds on average. 
Hydrogen bonds of the chains with water are obviously not counted employing an implicit 
peptide-solvent model and averages shown in the figure reflect only intrinsic hydrogen bonds. 
As can be seen in Fig. |9](a), the peptide with the randomly permuted sequence S3 behaves 
noticeably different than sequences SI and S2. For S3, (nn-bond) decreases more smoothly 
with increasing temperature, i.e., the breaking of the individual hydrogen bonds is a process 
of relatively weak cooperativity. In this case, the hydrogen bonds are comparatively weak 
and the two independent helical segments and, more globally, the helical phases, are not 
very stable, as the fluctuations show in Fig. Mjo). This effect is mainly due to the position 
of the proline in the chain. Although its distance from the ends is identical in all three 
sequences, the asymmetry with respect to the dihedral constraints destabilizes the larger 
helical segment of S3 (whose average length is smaller than the helices of SI and S2) and 
leads to a noticeable probability of forming a second small and weak helical segment. 

Finally, we have also calculated the radius of gyration as a rather global geometrical quan- 
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(a) 




FIG. 8: Probability profiles of each residue to be a part of a 3-residue helical segment as a function 
of temperature for (a): SI, (b): S2, and (c): S3 in solvent. 

tity, which is mainly useful for quantifying the structural collapse caused by a conformational 
transition. In contrast to the © collapse transition of polymers between non-structured 
globular and random-coil conformations, a helix-coil transition is rather a crossover from 
non-structured conformations to conformations with highly ordered segments (helices). For 
this reason, the gyration radius is a too rough measure for the order in the helical phase and 
is therefore of less importance for the understanding of secondary-structure formation and 
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FIG. 9: Average number of hydrogen bonds (a) and derivative (b) vs. temperature for the three 
peptides in solvent. 
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cannot be deduced from the short-range interactions 

In Figs.fTOlfa) and (b), the average radius of gyration and its fluctuations are shown. Being 
a measure of the compactness of the molecule, small values of the gyration radius indicate 
more tight-packed structures. For all sequences, the average radius of gyration changes 
from 6.0-6.5 A to 9.0-9.5 A. The peaks in the derivative of i?Gy indicate slightly higher 
conformational transition temperatures than the identified temperatures from fluctuations 
of energy and helicity. Sequence S3 possesses the most compact conformations in the helical 
phase - although the average number of hydrogen bonds is the smallest, as can also been 
seen in Table [Tll where the differences between maximal and minimal values of average 
radius of gyration and mean number of hydrogen bonds are given. 

In Table IIII[ we have listed the peak temperatures identified from fluctuations of several 
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FIG. 10: Average radius of gyration (a) and derivative (b) vs. temperature for the systems in 
solvent. 

thermodynamic quantities discussed in this section for the three peptides SI, S2, and S3. 
The peak temperatures of the peptides in solvent are compared with the corresponding 
transition temperatures identified for the systems in vacuum. While for the peptides in 
solvent the transition temperatures are relatively independent of the fluctuations considered, 
the deviations for the vacuum systems are noticeable. This is not surprising, as the systems 
in solvent are stabilized by the environment and the hclix-coil transition is a cooperative 
effect that is accompanied by a strong coupling to the solvent. In vacuum, the conformational 
freedom is much larger, and the helix-coil transition rather an entropic effect. In this case, 
the finiteness of the systems is more influential than for the peptides in solvent. 
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TABLE II: Differences between maximal and minimal values of average gyration radius A(i?Gy) = 
(i?Gy)max - (-RGy)min and mean number of hydrogen bonds A(nH_bond) = ("-H —bond) max 
('^H-bond)min for each Sequence over the whole temperature range. 



—bond/ 



SI 


2.6 


5.4 
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S3 
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4.1 



TABLE III: Helix-coil transition temperatures Tc in K, read off for the three sequences SI, S2, 
and S3 in vacuum and solvent from the various quantities discussed in the paper. The errors are 
in the range ±10 iT. 
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S2 Vac. 
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S3 Vac. 
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410 
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V. MULTICANONICAL HISTOGRAMS AND FREE-ENERGY LANDSCAPES 

For the study of the hehx-coil folding channels, it is useful to investigate the two- 
dimensional multicanonical histogram of the energy E and a suitable system parameter, 




FIG. 11: (Color online) Lowest-energy reference conformations of sequence S2 (left) and S3 (right) 
in solvent. 
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FIG. 12: (Color online) Multicanonical histograms of overlap parameter Q and energy E for the 
peptides (a) SI, (b) S2, and (c) S3 in solvent. 

which is chosen here to be the angular overlap parameter Q, as defined in Eq. (fTTl) . As 
reference conformations, required for the calculation of Q, we use the lowest-energy con- 
formations found in the multicanonical simulations. These structures are shown in Fig. [2] 
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(right) for SI and in Fig. [TT]for S2 and S3. 

The multicanonical histogram is obtained from the multicanonical time series 



(12) 



where the sum runs over the Monte Carlo steps t. The summation over Q yields the "fiat" 
multicanonical energy distribution. Since the conformational energy E, if replaced by the 
average energy {E), is directly related to the temperature, the histogram H{E,Q) contains 
sufficient information for the description of the simple folding transition. In Fig. [T21 we 
have plotted the multicanonical histograms for the three sequences SI, S2, and S3. In all 
three cases, we find a noticeable turn in the distribution from small values of Q, which 
correspond to random-coil conformations, to values closer to unity, where the conformations 
are similar to the helical reference conformation. In correspondence to the interpretation in 
the previous section, the transition is stronger for the wild- type and mutant sequences SI and 

52, respectively, and rather a two-step process in the case of the randomly permuted sequence 

53. The remarkable bifurcations at very low energies are possibly indications for metastable 
conformations which can be considered as weakly disturbed reference conformations. It 
should be noted that the angular overlap parameter is calculated by comparing all dihedral 
angles. This means that deviations in side-chain dihedral angles lead to Q values different 
from unity, although the backbone dihedral angles could take almost the same values. Helical 
structures are mainly due to cooperativity along the backbone, but differences in the precise 
side-chain locations do not destabilize the helical structures. Therefore, these kinds of glassy 
transitions, which happen under extreme conditions (very low temperatures!) are interesting 
but not in the main focus of this work. 

In Figs. [T3ra)-(c). we have plotted the free-energy landscapes for the three peptides. Here, 
we assume that the angular overlap parameter Q is a suitable measure for the structural order 
of the peptides. The free energy as a function of this "order" parameter and temperature is 
then given by 



F{Q,T) 



RT\np{Q,T), 



(13) 



with the distribution of the overlap parameters 




(14) 
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(a) 




FIG. 13: (Color online) Free-energy landscapes F{Q,T) for the peptides (a) SI, (b) S2, and (c) S3 
in solvent. 

where the integral runs over all possible conformations The free-energy plots in 
Figs. [T3r a)-(c) confirm that all three peptides experience a conformational transition follow- 
ing a single main folding channel. Although the helix-coil transition separates the pseudo- 
phases of random conformations and the long-range ordered helical phase, no noticeable 
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signal of pseudo-phase coexistence is observed, i.e., the transition appears rather second- 
order- than first-order- hke. 



VI. SUMMARY 



We have analysed thermodynamic properties and folding channels in the free-energy land- 
scape for three synthetic 12-residue peptides which exhibit remarkable adsorption affinities to 
semiconductors 5|, l6|. Employing an all- atom model based on the ECEPP/3 force field 25 1 
with OONS implicit solvation parameter set [21] and applying the implementation of the 



multicanonical Monte Carlo simulation method in the SMMP package [27|], we found in all 
three cases indications for a strong helix-coil transition. Independent of the fiuctuations 
studied, the peak temperatures are very close to each other - despite of the smallness of 
the peptides. Since experimental verification and biochemical structural analysis of these 
peptides are still pending, a comparison with experimental data is not yet possible. 

Our predictions for the transition temperatures are probably too high, as is to be ex- 
pected by using implicit-solvent models. Therefore, we expect that the helix-coil transitions 
could happen under reasonable environmental conditions such that it should be possible 
to verify our predictions experimentally. This is an important issue, since it is generally 
expected that selective synthetic peptides and polymers may play an essential role in future 
nanotechnological applications. 
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